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I Introduction 



The semileptonic process B ulu is known to provide information on the B ^ n 
form factor which for = is related to the factorized color-favored 

("T", or "tree") contribution to B^ vr+vr". In the present paper we show that 
while present semileptonic data have not yet reached adequate precision, with around 
500 B Triu events in which rate and spectrum are measured one can improve 
the accuracy of T by a significant amount. We then discuss the benefits of such a 
determination. 

A connection between the decays B^ — ^ 7r~£"'"z// and B^ tt+tt" was noted some 
time ago by Voloshin who derived the relation 



r(fiO ^ TT+TT-) 127r2/2 



13.7 (/^ = 131 MeV) 



using a pole model for the B ^ it form factor F^{q^). This relation assumes the 
dominance of a "tree" (T) contribution to B^ -k^-k^ in the notation of Ref. 
The CLEO and Belle 0] Collaborations have measured the branching ratio for the 
semileptonic process. Averaging their results yields 

^ 7r"e+z/e) = (1.4±0.3) X 10"^ (2) 

while an average of CLEO [Q, Belle [Q, and BaBar 0] (5° and i?°-averaged) branch- 
ing ratios ^ gives 

^ TT+vr-) = (4.4±0.9) X 10"^ . (3) 

The experimental ratio of these two branching ratios is T{B^ -K^e^Ve) /'^{B^ — ^ 
vr+TT^) = 32 ± 9, a factor of 2.3 above Eq. (|l|), which indicates either that the "tree" 
contribution is substantially overestimated in (|I]), or that some other process is inter- 
fering destructively with the tree amplitude to reduce the B^ —>■ 7r"*"7r" decay rate. A 
prime candidate for this amplitude is the "penguin" , or P amplitude in the notation of 
0. If this amplitude were sufficiently important to reduce the expected B^ —>■ 7i~^7i~ 
rate by roughly a factor of 2.3, it could have important effects on the extraction of 
the weak phase a = 02 entering the Cabibbo-Kobayashi-Maskawa (CKM) matrix [^]. 
This question has now acquired particular urgency as a result of the first report of 
results on CP-violating parameters in B^ — > vr+vr" . 

Many attempts have been made to use data to estimate the "penguin pollution" of 
the B^ TT+TT^ amplitude, including an isospin analysis requiring the measurement 
of B^ 71^71^ and B^ ti^tt^ decays |jTlI] (we assume charge-conjugate processes 
are measured when required), methods which use only a partial subset of the above 
information |]12|, and numerous methods based on fiavor SU(3) [|], |13[. Earlier data 



hinted that the penguin amplitude was interfering destructively with the tree in B^ 



14, 15 



In the present paper we describe measurements of -B'' ^ tt e+i^g decays which 
can significantly improve information on the magnitude of the tree (T) contribution 
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to — > TT+vr^. Such an improvement is needed to tell whether tree and penguin 
amplitudes are really interfering destructively in tt+tt". We discuss the role 

of the B* pole in this process, whose contribution is related through heavy quark 
symmetry to a recent CLEO measurement of the D*Dti coupling constant [1^. We 



then show how information on T helps to determine the weak phase a using limits 
on CP violation in B^ —>■ 7r"^7r~. 

Our approach differs from that advocated in Refs. 0, 0, |l8l, in which the tree 
amplitude is estimated from the rate for B^ — > tt+tt'^. In that process, there is an 
additional color-suppressed amplitude (called C in the language of Ref. whose 
magnitude and phase with respect to T cannot be independently estimated using 
present data but must be calculated. One then has A(B~^ tt+tt") = —{T + C)/^/2, 
and with C ~ O.IT, one arrives at estimates rather similar to those in the present 
paper. (The C amplitude was neglected altogether in Ref. [0].) The semileptonic 
process avoids dependence on the theoretical calculation of C/T. 

In Section II we give some basic expressions for the B^ Ti~e^Ve and B'^ ti^ti~ 
decays. Information on the B ^ ir form factors is reviewed in Section III. The 
D*Dtt measurement and its implications for the B^Bn coupling and the B* pole in 
the B ^ TT form factor are described in Section IV. We then bracket the possible 
magnitudes of the tree amplitude T depending on measurements of the spectrum in 
B^ — s> 7r~e"'"i/e (Section V). The extraction of the penguin amplitude from B Ktc 
decays with the help of flavor SU(3) allows us to determine the extent to which 
P and T are interfering destructively in B^ vr+vr", and hence to determine the 
correction to the weak phase a which is needed when extracting it from CP- violating 
asymmetries in that process (Section VI). We summarize in Section VII. 



II Semileptonic and nonleptonic tree decays 

For a generic heavy-to-light decay if — >■ vr, the non-perturbative matrix element is 
parametrized by two independent form factors: 

{7i{p)\u^,b\Hip + q)) =(2p + q- g ""^ ~ ) F+(g^) + g, ""^ ^ Fo(g^) , (4) 

with H being a. B or D pseudoscalar meson. The subscript H has been suppressed in 
the two form factors. In the case of massless leptons (which is an excellent approxi- 
mation for £ = e,fi), only F^{q^) contributes to the differential decay rate 

riv n"^ IT/ |2 

|^(if° .~£^v,) = "^j^lpW'lF^iq'r , (5) 

where VqQ is the relevant CKM matrix element. We will take |Vcrf| = 0.224 ± 0.016 
and \Vub\ = 0.0036 ± 0.0010 from Ref. [||]. To obtain the total width, one should 
integrate Eq. (Bl) over the entire physical region, 

< g2 < 

{tuh — rriT^) , which requires 

the precise knowledge of the normalization [i.e., -F+(0)] and dependence of the form 
factor. 
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Figure 1: Feynman diagrams for semileptonic and nonleptonic tree decays of a 
meson. 



The lepton pair can be replaced with a pion, as shown in Fig. |I| for the decay of a 
meson. The resulted diagram is the "tree" contribution to the nonleptonic decay 
B^ TT+TT^. In the limit of small m.,^, the two diagrams in Fig. |l] are related by the 



Bjorken relation pO | 



(6) 



where |ai| is the QCD correction. We shall take |ai| = 1.0, which is a sufficiently 
good approximation for our present purpose. 



Ill H ^ 7T form factors 

In the absence of a spectrum measurement, one cannot directly employ Eq. (^) to cal- 
culate T. Present extraction of T using this relation relies on assumptions of particular 
form factor shapes. One can test such assumptions using data on the B*Bt{ coupling 
extracted using heavy quark symmetry from the corresponding D*Dit coupling, and 
using present information from lattice gauge theories. Form factors parametrized in 
a manner consistent with such constraints can then be used to anticipate the number 
of events necessary to extract T from (^) in a model-independent way. 

Lacking experimental measurements of the form factors and -Fo(g^), people 

have proposed several models to describe their behavior, among which is the 



single-pole model: 



2\ fn* gH*Hn 



IrriH* 1 — q /mjj* 



where we adopt the following convention: 

{0\V,\H*{p,e)) = fH*mH*e,, (8) 
{H-{p)n+{q)\H*%p + q,e)) = Qh^-hM-^)- (9) 

However, this form factor gives total widths of n^i^Ui and B^ n^i^ug which 

are both larger than the experimental values, as will be shown in Section IV. So the 
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Table I: Vector meson decay constants (MeV) from different calculations. 



D* 



B* 



245 ± 201^ 
268tl^ 
327 ± 13 
354 ± 90 



196 ± 24+f 
236ti 
252 ± 10 
206 ± 39 
190 ±30 



Becirevic et al. 
UKQCD PI 
Hwang & Kim |7| 



Wang & Wu p 
Huang & Luo 



monopole form factors are not enough to describe the physics involved in the H ^ n 
decays. 

Multipole form factors naturally become our next choice. On the basis of lattice 
gauge theory calculations, Becirevic and Kaidalov proposed a simple parametriza- 
tion which is essentially a dipole for 



(1 - q^/mjj,){l - aHq^/mjj.) 
fotf) = , '"'-'-""l , ■ (11) 

In the infinite quark mass limit, the quantities {cHy/rnn, (1 — aH)mH, (Ph — 1)"^h) 
should scale as constants, ch is related to the coupling constant gH*H-K as 

Ch = — • 12 

This parametrization has enough freedom to describe lattice results, which typically 
are obtained for values of above about 13 GeV^ p2|, |2^. We shall employ 



it to judge the statistical accuracy needed in extrapolating the B rriu spectrum 
to = m^, where the Bjorken factorization relation (^) provides an estimate of T. 
A similar problem arises when one wishes to extrapolate to the zero-recoil limit in 
estimating the CKM matrix element \Vcb\ from the exclusive process B — >■ D^*Hu, 
since both the normalization and shape of the spectrum have to be determined. 

It should be pointed out that fn*, fs* and gB'B-K are far from being determined, 
though go'D-K has been measured [0. Very different values of fjj* and fs* have been 
obtained on the lattice and in various models (see Table | [^, ^ We 

will discuss gB*Bn in Section IV. 



IV Implications of Qd^d-k measurement 

We now describe the CLEO measurement of the D*D7r coupling constant and 



review its significance in the light of earlier predictions [^, 0, |3^. The observed 
value of the total D*^ width is T{D*^) = (96 ±4 ±22) keV, in satisfactory agreement 
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Table II: Predictions for decays D* Dtt and D* D^y based on comparison with 
K* Kit and K* K^y decays. 





Predicted 




Experimental 




Partial Width Branching 


Ratio 


Branching Ratio 


Decay 


(keV) 


(%) 


(%) 




25.9 


30.9 


30.7 ±0.5 




56.9 


67.8 


67.7 ±0.5 




1.1 


1.3 


1.6 ±0.4 




83.9 








39.7 


70.6 


61.9 ±2.9 




16.5 


29.4 


38.1 ± 2.9 




56.2 







with a prediction of 84 keV made some time ago by comparison with K* Ktt and 
K* — > Kj decays |3^. Other predictions of |0 are compared with the current ex- 



perimental situation in Table |1[ The agreement is not bad, and can be improved 
by assuming about a 30% increase in the absolute square of the matrix element for 
the magnetic dipole transitions D* —>■ Dj with respect to the value in Refs. The 
experimental branching ratios at the time of these predictions differed from them 
much more significantly. 

A more detailed set of calculations was performed on the basis of chiral and 
heavy quark symmetry [31|, taking into account SU(3) violating contributions of 



order m^^. The experimental values are consistent with the predicted correlation 
between B{D*^ D^l) and r(D*+), as shown in Fig. ||. 

The observed D*^ width can be related to a dimensionless 0*0-71 coupling con- 



stant g by the expression [31, 33 



r(D*+^D°7r+) = -^b-;|3 , (13) 

where /vr = 131 MeV and Ijvl = 39 MeV/c. Using the branching ratio in Table || we 
find T{D*+ D^TT+) = 65 ± 15 keV and g = 0.59 ± 0.07. Therefore 

gD'D. = ^g = 17.8 ±2.1 , (14) 

where m|) = 1973 MeV is the spin-averaged mass of the D^*^ meson. Taking this value 
of go* Dtt and fo* = 200 MeV (which is more than Icr smaller than any determination 
in Table |), we get B{D° vr-e+z/e) = (4.9 ± 1.2) x lO'^, still larger than the 
experimental value (3.7 ± 0.6) x 10^^ [0. Higher values of /d* yield even larger 
branching ratios. 

Heavy quark symmetry (HQS) predicts 

^B-B. = -7-^^ = 47.9 ± 5.7 . (15) 

In 
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Figure 2: Prediction of Ref. |^ for r{D*^) as a function of the branching ratio for 
D*~^ — > D^'j, including leading SU(3)-breaking effects. Lines show predicted bounds. 
The plotted point shows current data |TB|, . 



Again, even if we take a comparatively small value of fs* (=160 MeV) and assume 
a large (e.g., 40%) violation of HQS (so that gB*Bn can be as small as 29.0), we will 
get a branching ratio B{B^ 7r~e^z/e) = (2.6 ± 1.4) x 10~^ which is still larger than 
Eq. (^. Thus we are justified to suspect the single pole form factor (|^. 



V Information on T from semileptonic decays 

The Bjorken relation (|^) establishes a useful connection between the semileptonic 
decays and the nonleptonic "tree" decays. Ideally, dT{B'^ n^i^Uij/dq'^ at = 
provides the "tree" contribution to the branching ratio for B'^ vr+vr" (aside from 
QCD corrections, which have been found to be a few percent in related processes). 
However, in practice one must measure the semileptonic decay spectrum over a range 
of in order to accumulate a sufficient number of events, and therefore must model 
the spectrum shape, as in extracting \Vcb\ from the spectrum for B D'^*Hh'. 

The dipole form factor has enough parameters to allow modeling both a normaliza- 
tion and a spectrum shape. We use it to gain an idea of the statistical requirements 
for a useful spectrum measurement. The experimental branching ratio (0) for the 
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0.2 0.3 0.4 0.5 0.6 

Figure 3: The dependence of cb on as for given values of B{B'^ n^e'^Ue). Solid line: 
B{B° TT-e+z/e) = 1.4 X 10"^; upper dashed line: B{B° n'e+v^) = 1.7 x 10^^; 
lower dashed line: 7r~e+i/e) = 1-1 x 10""^. 



semileptonic decay — > 7r~e~^i'e puts a strong constraint on the dipole parameters 
cb and a^, as shown in Fig. Accordingly, the "tree" branching ratio for B^ —>■ tt+tt" 
is constrained to lie in a certain range (Fig. ^). It should be noted that Fig. | does not 
depend on \Vub\, though Fig. ^ can be altered by any change in \Vub\- We can always 
combine \Vub\ with cb and view |Vu6|cb as a single parameter. This observation plays 
an important role in estimating T from Fig. ^. 

To determine ub and hence cb and i3tree(-B° — * tt+tt^), one can measure the 
normalized spectrum (f ^) for 5° —> 7i~£~^h'i. Note that f ^ is independent of 
and \Vub\- Thus measuring its dependence on will give us very clean information 
about as- Fig. ^ shows us that a comparison of the spectrum in the interval < 

< 11 GeV^ with that for 11 < < 26 GeV^ should be useful in determining ob- 

A recent lattice calculation |^3| obtains values of ranging from about 0.2 to 
0.6, Cb from about 0.3 to 0.6, and -F+(0) around 0.27 with a 25% error. A more recent 
analysis from QCD sum rules on the light-cone obtains -F+(0) = 0.26 ± 0.08, in 
good agreement with the lattice result. This implies that parameters are within the 
ranges quoted in Figs. |^-|^, and leads to values of BtrcciB^ tt+tt^) ranging between 
about 4.5 X 10"*^ and 11 x 10"*^, as in Fig. |. 

Given the central value of B{B — > T^iu), Fig. |^ implies that an error Aas = 0.1 
will correspond to an error in AStreel-B" — > tt+tt") of about 10%, or an error in T of 
about 5%. An additional error will be associated with the statistical error associated 
with B{B — > Triu) itself. We shall determine the number of events required to achieve 
an error of Aa^ = 0.1, and estimate the corresponding total error in T. 



8 



1.0-10 



8.0-10 



o 
PQ 

(D 



h 6.0-10"^ 



PQ 



4.010 




0.6 



Figure 4: The dependence of i3tree(-S° — > tt+tt ) on as for given values of B{B^ 
TT-e+i/e)- Solid line: B{B^ T^'e+v^) = 1-4 x lO"^; upper dashed line: B{B^ 
TT-e+i/e) = 1-7 X 10-^ lower dashed line: ^ 7r-e+i/e) = 1-1 x IQ-^. 
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Figure 5: Normalized spectrum of 5° — * 7r~£~^i'£ for various values of as- At low 
the curves correspond to as — 0.20, 0.30, 0.40, 0.50, 0.60, from top to bottom. 
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Table III: Dependence of the fraction f of —>■ n i'^i'e events below = 11 GeV^ 
on the parameter ub- 





0.2 


0.3 


0.4 


0.5 


0.6 


f 


0.618 


0.595 


0.568 


0.538 


0.503 



In Table |T| we show the fraction / of 5° ^ vr i^Ui events below = 11 GeV^ 



as a function of a^- In order to obtain an error of Aa^ = 0.1, one has to determine / 
to a precision of A/ = 0.023. For a total of N events in the spectrum, the error in / 
is A/ = \J f{l — f)/N, which is about 0.5/V^ for / near 0.5. Thus, one needs about 
(0.5/0.023)^ ~ 470 B — > niu events to achieve this accuracy. Such a sample will be 
associated with an error in the overall B niu rate of 1/V470 ~ 4.6%. When added 
in quadrature with the 10% error in i3tree(-B° tt+tt^) associated with the spectrum 
shape, this leads to an overall error of 11% in |Tp or 5.5% in T. One will need 
considerably more than 470/ B{B — >■ niu) ~ 3.4 x 10^ B decays to obtain a sample of 
this size, since the efficiency of reconstructing the semileptonic decay is small (e.g., 
slightly under 2% at Belle ^). The Belle Collaboration has reported a signal of 107 
events on the basis of 21.2 fb~^, but the background (148 events) is larger than the 
signal, and the branching ratio is dominated by systematic error. Thus a sample of 
about 4.4 times the present size would be the minimum needed to achieve the stated 
goal, with a larger sample required if background levels are to be reduced. 

VI Information on P and its interference with T 

We shall use present and anticipated information on T based on the methods described 
in the previous section, and flavor SU(3) to obtain information on P from the 
mainly-penguin process -B^ K^tt^. In this manner we shall end up with an estimate 
\P/T\ = 0.26 ± 0.08, to be compared to the value of 0.259 ± 0.043 ± 0.052 quoted by 



Beneke et al. on the basis of a theoretical calculation. (The inclusion of weak 



annihilation contributions in |3^ raises this value to 0.285 ±0.051 ±0.057.) Improved 
input data will potentially reduce the error on this ratio considerably, allowing for 
an estimate of direct CP-violating effects in B^ —>■ tt+tt^ with less recourse to theory. 
Furthermore, if |T| turns out to be incompatible with the experimental magnitude 
of the amplitude A{B^ —>■ tt+tt^) = — (T ± P), we shall obtain a constraint on the 
product cos a cos 5, where a is the CKM phase discussed previously and 6 is the 
relative strong phase between tree and penguin amplitudes. Our discussion will be 
an updated version of that presented in P^ . 

We shall quote all rates in units of {B^ branching ratio xlO^). Thus, the average 
(0) of B^ 7r+7r~ branching ratios implies 

+ |P|2-2|rP|cosacos(5 = 4.4±0.9 (16) 

in these units. With Btrcc{B^ tt+tt^) ranging from (4.5 to 11) xlO~^ we then 
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TT ' TT 







estimate |T| = 2.7 ±0.6. This is identical to the value obtained from 
with additional assumptions about the color-suppressed amplitude. 

The penguin amplitude can be estimated from K^-k^ . The average of 

CLEO PI, Belle [B], and BaBar M branching ratios M gives 



^ K^Ti+) = (17.2 ± 2.4) X 10" 



(17) 



leading to |P'p = (17.2 ± 2A){t°/t+), \P'\ = 4.02 ± 0.28, where we use the lifetime 
ratio tb+/tbo = 1.068 ± 0.016 Here P' refers to the strangeness-changing 6 — > s 
penguin amplitude, which is dominated by the CKM combination VtsV^l- 

We now estimate the strangeness-preserving b d amplitude by assuming it to 
be dominated by the CKM combination VtdV^l. This may induce some uncertainty if 
the lighter intermediate quarks also play a role . (A slightly different definition of 
P is used by p5| , p9| and avoids this problem.) We find 



\P/F\ 



Vis 



All 



p-ir/l ~ 0.22(0.80 ±0.15) , |P| = 0.71 ± 0.14 



where A, p, and "i] are parameters HO] describing the hierarchy of CKM matrix ele- 



ments. Combining these results, we find only that —0.1 < cos a cos 5 < 1, so that 
destructive interference is possible but not established. Reduced errors on |T| and 
\P\ will be needed for a more definitive conclusion. In particular, given the present 
central values, reduction of the error on |Tp to 11%, as achievable with 470 B niv 
events, would allow one to infer the presence of destructive interference at about the 
2.8cr level. 

With our present estimates of \P\ and |T|, we then find \P/T\ = 0.26±0.08. Errors 
on this quantity can be decreased by improving the measurements of the branching 
ratio for B 7i£u, by measuring its spectrum, and by reducing the error on \l—p—iri\, 



which we have taken to be greater than in some other determinations [41 



The presence of the P amplitude can affect the determination of the weak phase a 



using CP- violating asymmetries in B 



decays. The BaBar Collaboration |T0 



has recently reported the first results for this process. The decay distributions (/_) 
in an asymmetric e~^e~ collider at the T(4S') when the tagging particle (opposite to 
the one produced) is a B^ (B^) are given by [O 



f±{At) - e-'^^^^il ± S^^ sin(AmrfAt) t C^n cos(AmrfAt)] 



where 



and 



Here 



Stt 



2Im(A, 
1± |A^ 



1- lA. 



1± |A. 



At, 



-2j/3 



TT TT 



v4(5° 

a(b' 



TT TT 



TT TT 



A(50 ^ TT + TT-) 



(19) 
(20) 

(21) 
(22) 
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50 100 150 

a 



Figure 6: Relation between ctes as measured using 5"^^ = sin(2Q;efr) and the weak 
phase a for \P/T\ = 0.26 and 6 = (sohd curve). Dot-dashed curves correspond to 
±1(7 errors on \P/T\. The dotted hne corresponds to P = 0. The sohd and dashed 
hues correspond to the central and ±la values of S^^r recently reported by the BaBar 
Collaboration (allowing also for error in \P/T\). We show only the range associated 
with the region of CKM parameters consistent with other measurements. 



where 5t and 6p are strong phases of the tree and penguin amplitudes. To first order 
in \P/T\, using /3 + 7 = tt — a and defining 6 = 6p — 6t, we then have 



A. 



e^'" ( 1 + 2^ 



e sm a 



In the limit of small \P/T\ and vanishing final-state phase S, the S.,, 
sin(2a;efr), where 



(23) 



term is just 



Oes — a + 



sm a 



(24) 



A plot of this relation for \P/T\ = 0.26 ± 0.08 is shown in Fig. The BaBar 



Collaboration has reported S^^ = 0.03l[J:i ± 0.11 on the basis of 30.4 fb^^. The 
corresponding central and ±lcr values of a^s and a are shown as the solid and dashed 
lines on the figure. 



12 




Figure 7: Relation between 5"^,^ and Ct^tt for fixed values of a (solid curves) and 5 
(dashed curves). The values of a range in steps of 10° from 50° (right) to 100° (left); 
those of 5 range in steps of 15° from —45° (bottom) to 45° (top). Here \P/T\ = 0.26 
has been assumed. 



To first order in \P/T\, the Ct^t^ term may be written 

C^^ ^ 2|P/T|sin5sina . (25) 



The BaBar Collaboration's value ^ C^^ = -0.25^1]:^^ ±0.47 is consistent with zero, 



as one might expect for a small final-state phase 6. This measurement in the future 
will serve mainly to constrain S, given the limited range expected for \P/T\ and sin a. 
Such a constrained value of S will then be useful in interpreting the flavor-averaged 
branching ratio (^ in terms of the tree-penguin interference discussed previously. 
The combined measurements of the flavor-averaged tt+tt^ branching ratio and 

the coefficients Sj^t^ and C^rTr, when combined with independent determinations of |T| 
and \P\, should allow us to dispense with the assumptions that the final-state phase 
S is small and that the weak phase of P is dominated by the top quark in the loop. 

An example is shown in Fig. ^ of how S-^t^ and C^j^ measurements can be used 
to constrain a and 6. Values extracted from such plots can then be checked for 
consistency with Eq. ([T6|) to check our assumption that the phase and magnitude of 
P is dominated by the top quark. 
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VII Conclusions 



We have discussed rate and spectrum requirements in B ^ iriui decays needed to 
reduce errors in the tree-amphtude contribution T to 5° — tt+tt^. Better knowledge 
of T can be combined with an estimate of the penguin amphtude P to see if destructive 
tree-penguin interference is occurring in — > tt+tt", and to evaluate the correction 
to the time-dependent CP asymmetry parameters S'tttt and C^^. Present data lead to 
the estimate \P/T\ = 0.26 ± 0.08 but substantial improvement will be possible once 
the semileptonic rate and spectrum (particularly near = 0) are better measured. 
We have estimated that at least 470 B — >• niu events (about 4.4 times the present 
sample size at Belle) are needed to reduce the error on T to 5.5%. For a near 90° 
we predict — a — (15 ± 5)°. Destructive tree-penguin interference in 5° tt+tt^ 
could be significant if a were closer to the lower limit of about 56° allowed by the 
present analysis. The form factor F+(g^) measured in B ^ Triui also can be helpful 
in estimating the "wrong-sign" amplitude in B — > D*tt decays 021 . 
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